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Importance of far-infrared mapping in a spiral galaxy: AKARI 
observation of M 81 
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ABSTRACT 

The importance of the far-infrared (FIR) mapping is demonstrated for a face-on spiral galaxy, 
M81, by analyzing its imaging data at 65, 90, and 140 |xm taken by AKARI. Basic products 
are the dust temperature map, the dust optical depth map, and the colour-colour diagram. The 
main features are as follows, (i) The dust temperature derived from the total fluxes at 90 |_tm 
and 140 [J.m reflects the relatively low temperatures seen in the interarm and spiral arms ex- 
cluding the warm spots, rather than the high temperatures in warm spots and the centre. This 
indicates that the total FIR luminosity is dominated by the dust heated by the general interstel- 
lar radiation field, (ii) The galaxy is more extended at 140 |^m than at the other shorter wave- 
lengths, which reflects the radial dust temperature gradient, (iii) The dust optical depth derived 
from the FIR mapping is broadly consistent with that estimated from the FIR-to-ultraviolet 
luminosity ratio, (iv) The FIR colour-colour diagram is useful to identify a 'contamination' of 
warm dust. The existence of small-scale warm star-forming regions is supported in the bright 
spots along the spiral arms. This contamination also leads to an underestimate of dust optical 
depth (or dust column density). 

Key words: dust, extinction — galaxiesTSM — galaxies: individual: M 8 1 — galaxies: spiral 
— galaxies: structure — infrared: galaxies 



1 INTRODUCTION 

L> '. 

J_( ' The far-infrared (FIR) emission from gal axies is dominated b y 
, the thermal emission from dust grains (e.g. lHarper&Lowlll973h . 

5^ ■ and is known to be a good indicator of star form ation rate (e.g. 

d ' |Kennicuttlll998l :ll noue, Hirashita. & Kamavall2OO0l) . Since the ma- 
jor heating source of dust grains is usually stellar radiation, the FIR 
colour, which reflects the dust temperature, is one of the most di- 
rect indicators of the interstellar radiation field (ISRF). Moreover, 
the FIR emission is usually optically thin in the ISM, which means 
that we can trace all the radiation over the entire column. 

Spatially resolved images of galaxies are required to study the 
dust properties in different regions of a galaxy and the characteris- 
tics of individual star-forming regions. Also, since only integrated 
global properties are observable for distant star-forming galaxies, 
it is important to clarify the relation between the global and local 
characteristics. Specifically we should distinguish whether the ma- 
jor dust heating comes from individual star-form ing regions or gen- 
eral ISRlEI(e.g. Walt erbos & Greenawaltl 1 9961) to clarify the major 



origin of the dust heating source responsible for the FIR emission. 
However, because the spatial resolution of FIR observation is gen- 
erally poor, only galaxies with large angular size can be well re- 
solved. 

M 8 1 is particularly suitable for our purpose because (i) it is 
a nearby big galaxy with a face-on geometry; (ii) it can be com- 
pared with another well-known spiral galaxy, the Milky Way; and 
(iii) it is within a field of view of the scan observation of AKARI, 
which is useful to quantify the global properties and to subtract the 
background^ 

Which wavelengths in FIR are necessary? At long wave- 
lengths (> 90 urn), the dust emission is mainly contributed 
by relatively large grains in radiative equilibrium with the stel- 
lar radiation field. However, an emission excess in shorter wave- 
lengths (< 65 (im) is usually observed, which comes from very 
small dust grains suffering from stoc hastic heating, that is, heated 
transiently to a high tem perature dAannestad & Kenvonl 1 19791 : 
iDraine & Anderson! Il985l) . In order to estimate the large grain 
temperature, which can be used to indicate the intensity of the 



* E-mail: hirashita@asiaa.sinica.edu. tw 

1 The general ISRF is defined as the average radiation field one would mea- 
sure from a location somewhere in a galactic disc where the radiation is not 
dominated by any particular star or star cluster ( Walterbos & Gre enawaltl 




2 In this paper, we use the AKARI data. Within the AKARI archives, M 101 
also satisfies our req uirements, but its FIS data are taken as part of the FIS 
calibration program ISuzuki et ai]|2007t) and are not open to public (i.e., a 
special treatment for the data is necessary). Other spiral galaxies observed 
by AKARI are too near or too far for our aim. 
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ambient stellar radiation, at least two bands are needed at wave- 
lengths > 90 \im. Furthermore, we require three wavelengths in 
the FIR to make a FIR colour-colour diagram. It has been shown 
that the FIR colour-colour relation of galaxies i s useful to derive 
the wavelength dependence of dus t emissivity JPaleetalJboOll ; 
iNagata et al.ll2002l : Iffibi et alj|2006h . The FIR colour-colour rela- 



tions in the Galactic disc and the Large/Small Magellanic Clouds 
are tight, indicating that the wavelength dependence of dust FIR 
emission can be considered to be uniform within and among those 
galaxies JShibai. Okumura. & Onakall 19991 ; Iffibi et alfcOOfj) . Such 
a common wavelength dependence of dust emission can be nicely 
repr oduced by a common grain optic al properties and size distribu- 
tion tairashita. Hibi. & Shibail2 007 ) . In this paper, the FIR colour- 
colour diagram also proves to be important for investigating the 
dust emission characteristics. 

Slow-scan data of Far Infrared Surveyor (FIS) onboard AKARI 
( TMurakami et alj|2007l : iKawada et ai1l2007l) with the 65, 90, and 
140 u.m bands is enough for the above requirements. (We do not use 
the 160 [im data because the data quality is not good enough for our 
purpose.) The 90 and 140 [im bands are used to estimate the large 
grain equilibrium temperature, and a colour-colour diagram can 
be made with the 65, 90, and 140 \im bands. Besides, by using the 
AKARI FIS data, a direct comparison amon g M 8 1 and a large num- 
ber of galaxies in the AKARI All-Sky data dYamamura et alj|2009T) 
can be made because of the common observational facility. Very re- 
cently0 S uzuki et al.l j2010h has reported an analysis of the AKARI 
data of M 81. They focus on the star formation law but do not in- 
vestigate the details of the dust temperature distribution within the 
galaxy and the FIR colour-col our diagram. M81 is also o bserved 
by the Spitzer Space Telescope dPerez-Gonzalez et alj20061) . whose 
70 \vca and 160 [im data can be compared with our data. How- 
ever, the dust temperature from the 70 \im and 160 u.m bands can 
overestimate the equilibrium temperature because of the possible 
contamination in the 70 [im emission by the stochastically heated 
very small grains. Thus, another FIR band at > 90 yim i s des irable. 
This galaxy has recently been observed by Bendo et"ai1 d2010h with 
Herschel Space Observatory jPilbrattetal]|2010h . whose long- 
wavelength coverage is particularly useful to examine the robust- 
ness of our analysis against the inclusion of long-wavelength data. 
Although the Herschel data are better than our AKARI data both in 
the spatial resolution and sensitivity, we concentrate on the AKARI 
data in this paper because not only the data of M 8 1 is available but 
also the All-Sky Survey data can be used to test if the FIR colours 
obtained for M 81 can be discussed in terms of a large number of 
other nearby galaxies taken by the same facility. 

In this paper, we aim to point out the importance of FIR map- 
ping data by deriving some fundamental quantities about dust (dust 
temperature, dust optical depth, etc.). At the same time, we try 
to clarify the relation between global (integrated) quantities and 
spatially resolved quantities. Some important features in the FIR 
colour-colour diagram of a spatially resolved image is also focused 
on. The method used in this paper is general and can be applied 
to future data taken by Herschel or the Atacama Large Millime- 
tre/submillimetre Array (ALMA)Q 

This paper is organized as follows. We explain the data anal- 
ysis in Section|2] and describe some basic results related to global 
properties such as morphology and radial dependence in Section 
[3] We investigate the distribution of dust temperatures and the 



colour-colour diagram in Section [4] We discuss the results ob- 
tained for M 8 1 in general contexts in Section [5] Finally, Section 
[6] gives t he conclusion. The di stance to M81 is assumed to be 
3.63 Mpc dFreedman et all 19941) throughout this paper. At this dis- 
tance, 1 arcmin corresponds to 1.06 kpc. 



2 DATA ANALYSIS 

M81 was observed by AKARI/F1S as a pointed observation (PI: 
FIS Team). The data are taken from Data Archives and Transmis- 
sion System[f| Photometric observations were performed with four 
bands: N60 (central wavelength: 65 iun), WIDE-S (90 Jim), WIDE- 
L (140 tun), and /V160 (160 (J.m) with an observational mode of 
FIS01 (photometry/mapping mode), a reset interval of 0.5 s and a 
scan speed of 15 arcsec s _1 . We utilize the images of the 65 |i.m, 
90 (im, and 140 [im bands. The quality of the 160 u.m band image 
is not good enough for our purpose. The FIS observation of M81 
is composed of two pointings: one covers almost all the M 81 area 
except for the south-east edge, which is observed by the second ob- 
servation. We combine those two images to obtain the entire M 81 
image after the background subtraction (Section[2j3j, and the anal- 
ysis procedures before combining the images are identical for the 
two images. The m easured FWHMs o f point spread function are 
37", 39", and 58" dKawada et al]|2007l) . which correspond to 650 
pc, 690 pc, and 1000 pc, respectively, at M 81. Thus, internal struc- 
tures such as spiral arms can be clearly identified on the images 
(Sectionl3~lY 



2.1 FIS Slow-Scan Tool 

The raw data we re reduced by using the FIS Slow-Sca n Tool (ver- 
sion 20070914: IVerdugo. Yamamura. & Pearson 2007)0 The pro- 
cess includes flagging of bad data, measurement of sky signal, 
dark and response correction, flat-fielding, and construction of co- 
added images. We used the local-flat for the flat-fielding. The out- 
put image grid size is chosen to be 30", which is about half of the 
beam size of 140 u.m. We confirmed that the following results are 
not sensitive to the selection of the grid size. The final image is 
smoothed with a length of ~ 40" for the SW (65 [im and 90 u.m) 
images and ~ 60" for the LW (140 u.m) image to avoid the effect 
of small scale detector noise. 

It is known that the AKARIfFIS detectors un derestimate the 
total flux probably because of the slow response IShirahata et al.l 
2008). Thus, we multiply the correction factors for the intensity, 1 .7 
for 65 [im and 90 [im, and 1.9 for 140 u.m, are multiplied to each 
band. Because the factors are similar to all the bands, the colours 
(flux ratios) and dust temperatures are not significantly affected by 
this correction. 



2.2 Colour correction 

Since the intensity at the central wavelength is derived by assum- 
ing a flat spectrum (vl v — constant) in the FIS Slow-Scan Tool, 
a colour correction should be applied. We applied correction to the 
WIDE-S (90 tun) and WIDE-L (140 (im) fluxes, assuming a spectral 
shape of I v oc B v (Td), where T<j is the dust temperature, /3 is 



3 After the submission of this paper. 

4 http://www.almaobservatory.org/ 



° http://darts.isas.jaxa.jp/ 

6 http://www.ir.isas.ac.jp/ASTRO-F/Observation/ 
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the emissivity index, and B v (T) is the Planck function. The emis- 

sivity index (3 is chosen to be 2 in this paper unless otherwise stated. I 1 1 1 1 1 

The value of ft is dependent on the grain compositio n, and = 2 is 

approp riate for astronomical silicate and graphite in lDraine & Led 69°09' - . S2 

dl984f) . We apply T d = 27 K and 19 K for WIDE-S and WIDE- 

L, respectively; the former (latter) temperature is derived from the t j ^ 

total fluxes of M81 at 7V60 and WIDE-S (WIDE-S and WIDE-L) ^. 69 06 - ^ 

(see Section [3T2T >. Consequently, the correction factor is 0.92 and o -it 

0.94 for the WIDE-S and WIDE-L fluxes, respectively, and the flux 3. 69°03' - 12 '^B^fc 

obtained by the Slow Scan Tool is divided by these factors. The 10 

uncertainty caused by the colour correction (5 per cent for WIDE-S 

and 1 per cent for WIDE-L for the temperature range observed in 69 00' - 

the M 81 disc; Section BTt is smaller than the errors caused by the ,S 1 

background fluctuation. For the narrow band N60 (65 iim) we do 68°57' ( a ) 

not apply the correction, since the colour correction only changes > 1 1 1 1 1 

the flux by less than 3 per cent. 09 h 56 m .5 56 m 55 ra .5 55 m 54 m .5 

a (J2000) 



2.3 Background subtraction 

Bright stripes along the scan direction caused by glitches or non- 
uniform detector sensitivity are observed in all the three bands. In 
order to eliminate these structures, the background levels are esti- 
mated for each line along the scan direction by averaging the in- 
tensities in the two sections before and after the scan of the M 8 1 
main body with a separation of about 18' and a total length of about 
20', and then subtracted. The RMSs of the background are 1.3, 0.7, 
2.9 MJy sr _1 for the 65, 90, and 140 urn bands, respectively. These 
values are estimated before eliminating the stripes to show the orig- 
inal uncertainty including the non-uniform sensitivity and the time 
variability of the detector response. 



2.4 Matching the positions 

We matched the positions of the images taken by two detectors, SW 
(N60 and WIDE-S) and LW (WIDE-L) according to the position of 
central peak of the galaxy image. The uncertainty in the relative 
position between the detectors is well below the grid size. 



2.5 The physical quantities derived from the data 

The intensity (surface brightness) at a wavelength A (frequency 
v = c/X, where c is the light speed) in each grid is denoted as 
I v (A). The intensity ratio, 7„(Ai)/I„(A2), is called Ai — A2 colour 
in this paper. We derive the dust temperature T c \ in each grid as 
follows. Since the dust emission is optically thin in the FIR (Sec- 
tion IQV Iu(X) = r(A)B„(Td), where r(A) is the optical depth 
at wavelength A. The functional form of the optical depth is as- 
sumed to be t(A) = Av 2 (y = c/X) with an unknown constant 
A, which is independent of the frequency. This functional form 
for the optical depth indicates that the frequency dependence of 
the dust emission coefficient follows oc v 1 with p = 20 Then, 
with two wavelengths (Ai and A2) selected, a set of equations, 
/„(Ai) = t(\x)B v (T a )\ m and 7„(A 2 ) = r(A 2 )B„(r d )U 2 , is 
solved to obtain Td and A. The dust temperature Td determined 
from wavelengths of 90 \im and 140 \im represents the temperature 



7 Although we consider a variation of f3 as a function of wavelength in 
Section 14.21 for a detailed theoretical model, the simple and common as- 
sumption of ft = 2 is useful to compare our observational results with the 
quantities obtained in the literature. 
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Figure 1. Images of M 8 1 . (a) The 90 urn image. The circles define various 
regions used for the analysis of the colour-colour relation in Section l4~2l the 
central region (C), the spiral arm regions (SI, S2, and S3), and the interarm 
regions (II, and 12). (b) The Ha image taken by Chen g et"aT] j 19971) (grey 
scale) overlaid with the 90 urn brightness contours. The highest contour 
shows 0.71 times the peak intensity and the intensity ratio of two adjacent 
contours is set to be 0.71. 



of large g rains in radiative e quilibrium with the ambient radiation 
field (e.g. lLi & Drainell200lh . and it is denoted as Tlg- 

It is convenient to convert r(A) = Av 2 to a commonly 
used indicator of dust optical depth. We choose Av (the extinc- 
tion in V band in units of magnitude) for such an indicator. We 
adopt the Galactic e xtinction curve for the conve rsion from t(A) 
to Av as derived bv lWeingartner & Draind J2001I) for Rv =3.1: 
A v = C{X)t(X), where the factor C(A) is 8.3 x 10 2 for 90 [im and 
1.9 x 10 3 for 140 irm. Although we derive Av by using the 90 [im 
values, we obtain the same value for Av also from the 140 [im val- 
ues because C(A) oc A 2 holds within 10 per cent between 90 and 
140 |^m. 



3 GLOBAL PROPERTIES 
3.1 Morphological features 

In Fig.[T] we show the 90 p.m image. The central bright peak, the 
spiral structures, and the spots along the spiral arms are clear. In or- 
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Table 1. Observational fluxes in units of Jy for the entire M 8 1 (global) 
and the individual regions shown in Fig.[T] 
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Figure 2. FIR SED of M 8 1 . The filled squares show our measurements, 
while the others are taken from the literature: The open diamon ds at 60 u.m 
and 100 u.m are taken from the IRAS data iRice et alj 1 19881). the open 
squares at 70 urn and 1 60 urn from the Spitzer data IPerez-Gonzalez et alj 
l2006t lDaleet alJl2007h . and the open triangles from the Herschel data 
teendo et alJl2010hT The solid (dotted) line shows the single-temperature 
SED which explains our measurements at 90 urn and 140 urn (65 p and 
90 urn). The dashed lines show the SEDs at the temperature of upper and 
lower bound of Tlq estimated by using the upper and lower fluxes at 90 urn 
and 140 urn. 



der to examine the correspondence between the FIR brightness and 
the star formation activity in M 8 1 , we compare our FIR image with 
the Ha image taken fro m the NASA/IPAC Extragalactic Database 
(NEDJ3 (originally from lCheng et alJI 19971) . The FIR emission re- 
ally traces the Ha central peak, the Ha spiral arms, the H II regions 
along the spiral arms. This confirms that the FIR emission is a good 
tracer of the star formation activities in galaxies. 



3.2 Spectral Energy Distribution 

The measured total fluxes, f„, at 65, 90, and 140 urn are shown in 
Fig.|2]and Table[T]together with the data taken from the literature. 
The fluxes are estimated by integrating intensities over a rectangu- 
lar region of width 15' and length 18' covering the galaxy main 
body. We conservatively adopt flux errors of 10 per cent for the 
65 ftm band and the 90 [im band and 20 per cent for the 140 (im 
band dHirashita et al.l R 



2008). As shown i n Table [T] we obtain the 



Suzuki et aljfco loh within the errors 



fluxes consistent with 

We also show two spectral energy distributions (SEDs) which 
fit the data at 90 and 140 vim and those at 65 and 90 |i.m with 
the functional form given in Section 1231 (i.e. /„ oc v 2 B V (T&)). 
We observe that the SED which fits the 90 \im and 140 \im 
fluxes underproduces the 65 urn flux. This indicates the contri- 
bution of very small grains (VSGs) to the emission at < 65 |im 
dAanne stad & K envonl 1 19791 : iDraine & Andersonl Il985h : Small 
grains with radii < 10 nm are stochastically (i.e. not in radiative 
equilibrium) heated to a high temperature, contributing to the in- 
tensity at the short (< 65 (im) wavelengths. The dust temperature 
Td derived from the 90 p.m and 140 p.m fluxes is 18.6 K, while that 



http://nedwww.ipac.caltech.edu 



Wavelength 


65 u.m 


90 urn 


140 urn 


Global 


94 ±9 


177 ± 18 


419 ± 84 


Global (S 10) a 


70 ±20 


189 ± 40 


356 ± 140 


C 


9.21 


12.8 


16.3 


SI 


3.61 


6.05 


10.7 


S2 


3.42 


6.52 


13.0 


S3 


1.65 


2.71 


5.32 


11 


2.06 


4.80 


10.8 


12 


2.12 


5.07 


11.7 



a Measurements by Suzu ki et all J2010h . 

Note. The conservative errors are 10 per cent for 65 and 90 urn, and 20 
per cent for 140 urn. 

derived from the 65 (im and 90 |J.m is 27.3 K. The former tempera- 
ture is interpreted to be the one of large grains which are in radiative 
equilibrium with the ambient ISRF (Tlg), while the latter does not 
present any physical dust temperature because of the contamination 
of stochastically heated very small grains to the 65 p.m flux. To es- 
timate the uncertainty in Tlg, we use the upper (lower) and lower 
(upper) values for the 90 and 140 p.m fluxes and obtain the upper 
(lower) temperature as shown in Fig. [2] 20.5 and 17.0 K for the up- 
per and lower temperatures, respectively. Thus, we finally obtain an 
estimate for the large grain temperature as Tlg — 18.6* \'g K from 
the global SED. We observe that the Spitzer and Herschel data at 
A > 160 p.m are broadly consistent with our SED fit to the 90 (im 
and 140 yim data within the error s. 

fpIrez-Gonzal ez et al. I d2006t) have shown that the Spitzer data 
can be fitted with two dust temperatures: their colder component 
with a temperature of 18 ± 2 K has a consistent temperature with 
our Tlg, while their warmer component with 53 ± 7 K shows a 
much higher temperature because they also used the 24 [im data. 
The temperature e stimated by the Herschel data is 17.8 ± 0.6 K 
jBendo et al]|2010h . Our estimate of the large grain temperature is 
consistent with the Spitzer and Herschel results within the errors. 

As a global quantity, we derive the total infrared (TIR) lumi- 
nosity (the integrated luminosity between 8 and 1000 p.m), which is 
used later in Section |4~T| For this aim, we first estimate the total flux 
traced by AKARI FIS . The AKARI FIR flux, /akari is defined by 
Hirashita etal] d2008l) to estimate the total flux in the wavelength 
range of 50-170 [im: 

f akari = fu(65 \im)Av(N60) + /„(90 \im)Au(WIDE-S) 
+/„(140 \i3n)Av(WIDE-L), (1) 

where / v {\) is the flux per unit frequency at wavelength A, and 
At; (band) denotes the frequen cy width covered by the band. Ac- 
cording to lKawada et"ai] d2007l) . we adopt At/(N60) = 1.58 THz, 
Av{WIDE-S) = 1.47 THz, and Au(WIDE-L) = 0.831 THz. With 
/akari, the AKARI FIR luminosity, L akari, is estimated as 



Lai 



= 4irD 2 f A 



(2) 



where D = 3.63 Mpc is the distance to M81. We obtain 
L akari = (3.1 ± 0.6) x 10 9 L . Considering that the 10-20 
per cent errors in the flux, we put an error of 20 per cent also to the 
luminosity. 

Next, we relate the AKARI luminosity and the TIR luminos- 
ity. There is a good correlation between L akari and Ltir (to- 
tal infrared luminosity emitted by dust) in the AKARI FIS Bright 
Source Catalogue, the first primary catalogue from the AKARI All- 
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Sky S urvey dYamamura et alj[2009l) . as shown by iTakeuchi et al.l 

boioh : 

log Ltir = 0.940 log Lakari + 0.914. (3) 

With this empirical formula, Ltir = (6.9 ± 1.4) x 10 9 Lq for 
M81. This luminosity is divided by 4%D 2 to obtain the total in- 
frared flux, Frm = 1-7 X 10 -8 erg cm -2 s _1 . Considering that 
the flux may be uncertain by 20 per cent, it would be appropriate 
to put a 20 per cent error for each flux or luminosity. Suz uki et al.l 
(2010) obtained smaller luminosity (4.0 ± 0.2) x 10 9 L for the 
total FIR luminosity. They consider contributions from two com- 
ponents: a cold component with a temperature of 22 K and a warm 
component with a temperature of 64 K, but these two components 
do not include the contribution from the emission at < 20 um. 

We can also check the consistency of the dust mass with the 
Herschel data. The total dust mass, Ma, is estimated as 



M d = 



f v (\)D 2 
K v B v (T d ) ' 



(4) 



where k v is the mass absorption coefficient at a frequency v. Us- 
ing the measured flux at 90 um, the dust temperature derived 
from the 90 and 140 um dat a, and k v = 34 cm 2 g _1 at 90 um 
dWeingartner & Draind l20oTh . we obtain M d = (3.2 ± 0.6) x 
10 7 Mq, where an error of 20 per cent is applied because of 
the flux uncertainty. This value is consistent with th e dust mass 
(3.4 x 10 7 Mq) obtained by the Herschel observation (Ben do et al.l 



3.3 Radial dependence 

In Fig. [3] we show the radial profiles in the three bands. The radial 
distances (r) are deprojected by considering the position angle (157 
degree) and the inclination angle (57 degree) taken from NED. The 
inclination angles given in NED range from 55 to 63 degree. There- 
fore the deprojected radius is uncertain by 10 per cent. The radius 
from the centre is divided into bins with a width of 50 arcsec, and 
the intensities of the pixels contained in each radius bin is averaged 
to obtain the intensity as a function of radius. The standard devia- 
tion for each radius bin is also shown with the bar. The flux errors 
are much smaller than the standard deviations. 

The radial profile can be divided into two components: bulge 
(r < 150 arcsec) and disc (r 150 arcsec). The 65 and 90 um 
band images show a bright bulge where the intensity increases 
rapidly as the radius decreases, while the 140 um image does not 
have such a sharp rise. The difference in the 'sharpness' of the bulge 
component among different bands is due to the high dust tempera- 
tures, which reflect an intense radiation from the active galactic nu- 
cleus (AGN) or from the stars with high surface density. The bulge 
component i s not prominent also in the long-wavelength Herschel 
SPIRE data dSauvage et alj20ld) . 

For the disk part, all the intensities in the three bands 
show a linear slope along the radius. The linear slope 
means that the dust emission is radially more extended than 
the stellar emission, whose radial pro file drops exponentially 
(e.g. | Baggett. Baggett. & Andersonll998h . Although lSauvage et~al1 
(2010) apply a fitting with exponential or Gaussian function to the 
FIR radial profile of the Herschel data, the fit is not necessarily 
good. The dispersion shown by the bars is relatively large around 
r ~ 400 arcsec. This radius corresponds to the radial extent of the 
spiral arms. 

To compare the radial extents in the three bands, we show 
the radial profiles of the colours, J„(140 um)//„ (90 um) and 




200 400 600 

Radial distance [arcsec] 

Figure 3. Intensity profiles in various wavelengths. The diamonds, trian- 
gles, and squares represent the averaged intensities as a function of depro- 
jected radial distance in the 65, 90, and 140 |j.m bands, respectively. The 
bars on the individual data points represent the standard deviation, not the 
eiTors. The solid lines are the fitting results with linear slopes for r > 150 
arcsec. 



T(65 urn)/ T(90 um) in Fig. [4] The former rises as the radius 
increases, showing that the 140 um image is more extended than 
90 urn, while the latter dose not show such a clear trend. The slight 
enhancement of the 65 um — 140 um colour and the plateau of the 
140 um — 90 urn colour around r = 6 arcsec are caused by the 
relatively high temperature in the spiral arms. 

Because of the linear-slope behaviour of the radial profile, the 
radial extent is well defined by the intercepts on the x axis in Fig. [3] 
The radial extents defined in this way are 706, 674, and 702 arcsec 
for the 65, 90, and 140 um images, respectively. The reason why 
the emission is more extended in 140 um than in 90 um can be 
the radial gradient of Tlg. The 140 um — 90 um colour is directly 
converted into Tlg by assuming a functional form of ^ 2 -B„(Tlg) 
for the SED. While 140 um is closer to the intensity maximum 
of the SED, 90 um at the Wien side is more sensitive to the tem- 
perature change, and decreases more rapidly outwards. A positive 
160 um — 70 um colour gradient and a negative temperature gradi- 
ent of cold dust component a re also shown with the Spitzer data by 
Perez-Gonzalez et al. ( 200^). 

iBendo et al.l d20ld) have also shown the radial profile of the 
70 um — 1 60 um colour by using the Herschel data. Their results do 
not show a monotonic decrease of the the 70 um — 160 um colour 
as a function of radius. The non-monotonic behaviour is similar to 
our AKARI 65 um — 90 um colour. On the other hand, the Herschel 
colours composed of the wavelengths longer than 160 um show a 
monotonic trend, which can be interpreted as a negative tempera- 
ture gradient. This is consistent with our 140 um — 90 um colour. 
Therefore, if a wavelength shorter than ~ 70 um is included, the 
colour behaves in a more complicated way possibly due to a rela- 
tively strong contribution from the local heating source to the short 
wavelength range. 



4 DETAILED FEATURES 

4.1 Distribution of dust temperature and optical depth 

The large-grain temperature (Tlg) is determined at each grid by 
the method described in Section 12.51 We only adopt the pixels 
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Figure 4. Radial profiles of the 140 \mi — 90 colour (upper) and the 
65 |rm — 90 |rm colour (lower). In the upper figure, the colour is converted 
into the equilibrium dust temperature Tlg by assuming oc i/ 2 B j/ (Tlq) 
for the SED shape. 



with brightnesses above 4 a of the background fluctuation, that is 
J„(90 nm) > 2.8 MJy sr" 1 and 7„(140 fim) > 11.6 MJy sr _1 
for the temperature map to avoid artificial features. The result is 
shown in Fig. [5] The spatial distribution of Tlg shows some struc- 
tures in M 81, especially the central peak and the warm spots along 
the spiral arms. There is an overall rough trend that the tempera- 
ture decreases outwards, which is consistent with the colour gra- 
dient shown in Section [331 The temperature is typically 25 K at 
the centre, 20-21 K at the warm spots along the arms, and 17-20 
K in the other diffuse regions, indicating more intense radiation 
field in the centre and in the spots. Moreover, the dust tempera- 
ture derived from the total fluxes (18.6 K; Section lT2l is consistent 
with that in the diffuse regions, which supports the view that dif- 
fuse (g^neraTj_rSRFjsjesr^sible for the major part of dust heat- 
ing |Walterbos & Greenawaltlll996h . The diffuse nature of the cool 
(Tlg — 18 K) dust is cons istent with the Spitzer 160 prn image 
dPerez-Gonzalez et alj2006h . 

We also present the spatial distribution of dust optical depth 
in Fig. [5] The optical depth at 90 [im is converted to the F-band 
extinction (Ay) (Section[23J. There is a weak anti-correlation be- 
tween Ay and Tlg- The large excess of Av in the north-west and 
south-west parts of the disc may be strongly affected by the con- 
tamination of stripes caused by unstable detector response. In Fig. 
[5] the most prominent stripe is shown by the dotted lines indicat- 
ing the scan direction. In the temperature map, there are two high- 
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Figure 5. Spatial distribution of (a) large-grain temperature (Tlg) an d 0>) 
dust extinction in the V band (Ay ). Both quantities are estimated from the 
90 )j.m and 140 [im intensities in the individual grids. The correspondence 
between the colours and the physical quantities is shown in the colour bar 
in each panel. The dotted lines indicate the most prominent stripe and the 
scan direction. 



temperature islands connected with a bridge on this stripe. The 
neighbouring region around this line with a width of 3 arcmin may 
not be reliable because of some stripes. Thus, future observations 
of the structures around the stripes are required. 

Our estimate of Ay is reasonable if the dust in each pixel has 
an uniform temperature. However, in the spiral arms, as shown in 
the Ha image (Fig.Q~|l, star-forming regions whose sizes are smaller 
than the grid size of the FIR image are clearly seen. If the dust 
temperature is biased to these regions, the dust optical depth can 
be underestimated. In fact, as shown later in Section |4~2l the FIR 
colour-colour relation also supports a significant contribution of 
warm star-forming regions to the FIR emission in spiral arms. The 
central part does not show a strong indication of such inhomoge- 
neous radiation field; that is, the dust temperature in that region 
might be well approximated by a single component. Thus, there 
is still a possibility that small Av in the central part is due to a 
deficiency of dust. However, since Ay is very sensitive to the de- 
termination of the dust temperature and our results are affected by 
unstable detector response, all the features in Ay should be con- 
firmed by future observations. 
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The Av value is roughly 0.5-2 in the M 81 disc. iBuat et all 

J2005I) show that FIR to ultraviolet (UV) luminosity ratio of galax- 
ies is related to the dust extinction because the efficiency of UV 
reprocessing into FIR by dust is related to the dust extinction. The 
UV flux of M 81 at the FUV band (1528 A) of the Galaxy Evolu- 
tion Explorer (GALEX) is F v = 179 mjy foale et alj|2007h . Thus, 
Ffuv = vF v — 3.51 x 10 -9 erg cm -2 s _1 . The extinction at the 
FUV band can be estimated b y using the fitting formula derived 
from the model calculations bv lBuat et alj 12005): 

A FUV = -0.0333y 3 + 0.3522?/ 2 + 1.1960y + 0.4967, (5) 

where y = log(FTm/ 'Ffuv)- Using the value for Ftir in Section 
13.21 we obtain Afuv = 1.46. This mat ches the mean value for the 
UV-selected sample in lBuat et al] j2005h in the local Universe. Now 
using the Galactic extinction curve given by Weingartn er & Draind 
feOOlh with R v = 3.1 (Av/A FU v - 0.389), we obtain A v = 
0.57. The extinction shown in Fig.[5]is the one over the entire col- 
umn in the galactic disc. If the stars are on average located in the 
middle of the disc thickness, the stellar extinction would be half of 
the values obtained in Fig.(5j that is, Av — 0.25-1. Ay — 0.57 is 
within this range. 



4.2 Colour-colour diagram 

The relation am ong the intensities at the three bands is investigated 
here. Following Ifflbi et al.1 ^2006), we investigate the FIR colour- 
colour relation. Fig. [6] shows the relation between the 65 [im — 
90 p.m colour and the 140 u.m — 90 u.m colour for the individual 
grids whose intensities are above 4 a of the background fluctuation 
for all the bands (i.e., the same criterion for the temperature map 
in Section |4~TV The overall trend from the lower right to the upper 
left can be interpreted as a sequence of dust temperature. 

Now we theoretically quantify the observed colou r -colou r 
relation by adopting the calculations in Hirashita e t al.1 j2007t) . 
who treat the dust heating by an ISRF and the dust cooling by 
thermal radiation to calculate the temp erature distribution func- 
tion using the framework developed by iDraine & Lil j200lh . The 
physical quantities that explain well the dust emission proper- 
ties in the solar neighbourhood a re adopted: the ISRF SED by 
Mathis. Mezger. & Panagial 19831), the grain size distribution by 



Mathis. Rum pl. & Nordsieck 1 19771) . and the heat capacity of grain 
materials bv lDraine & Lil ( feOOlh . The e mission coeffic i ent of dust 
(silicate and graphite) are taken from IDraine & Lee] |l984) for 
A ^ 100 urn and extrapolate d by assuming a functional form pro- 
posed by iReach et alj d 19951) (i.e. ft sm oothly changes from 1 to 
2 around A ~ 200 (J.m). As shown by iHirashita et alj J2007I) . a 
slight change in /3 affects the colour-colour sequence significantly. 
To check the consistency with other nearby galaxies in the colour- 
c olour diagram, adoptin g the same emission coefficient as adopted 
in IHirashita et all (2007 ) is crucial here. We vary the ISRF with the 
spectral shape fixed, and denote the ISRF i ntensity relative to the 
solar neighbourhood value as \. We refer to lHirashita et al.1 J2007I) 
for the details of the framework and some basic results. 

We show the FIR colour-colour relation for various ISRF in- 
tensity x m Fig- IU We only show the results for graphite, since 
silicate f ollows the almost iden tical FIR colour-colour relation to 
graphite (Hirashita et al. 2007). We observe that the FIR colours 
obtained for the individual grids are roughly explained with \ = 1— 
30, although most of the points are located systematically above the 
theoretical predictions on the diagram. The lower values of \ cor- 
respond to the lower- Tlg regions such as interarm regions, and the 
higher values to the centre and the bright spots in the spiral arms. 



The above theoretical colour-colour sequence is correct if the 
radiation field in a grid is approximated to be uniform. As clearly 
seen in the Ha image in Fig. Q] there are small-scale star-forming 
regions, which should host warmer dust because of high radiation 
field intensity. In order to examine the effect of such a 'contami- 
nation' of warm dust, we show the FIR colour-colour relation by 
mixing the results with x = 1 an d x = 100: the former value rep- 
resents the general ISRF, while the latter is taken as a representative 
high radiation field value. The fraction of the latter (i.e. higher x) 
component is denoted as fh\ that is, the intensity is calculated by 



h = (i- h)h{x = 1) + hh{x = 100), 



(6) 



where I v (x) is Iv as a function of x- The result is shown in Fig. [6] 
A slight contamination of the higher- x component with a fraction 
of /h = 0.003 significantly lift the colour sequence upwards, ex- 
plaining the upper part of the FIR colour-colour relation of the 
individual grids. This is because the 65 [im intensity responds 
most sensitively to the higher x component. The shift of the FIR 
colour-colour relation by the contaminatio n of a higher x. com- 
ponent is consistent with the conclusio n by iHibi et al] j2006l) and 
IHirashita et"ai] ( l2007l) . lHibi et al J J2006h called this upper sequence 
'sub-correlation', and the contamination effect 'overlap effect'. 

We also compare our results with the FIR colours of the galax- 
ies in the AKARI FIS Bright Source Cat alogue in Fig. [6] Th e anal- 
ysis of these galaxies has been done by IPoIIo. Rvbka. & Takeuchl 
(2010). Since there are an enormous number of galaxies, we only 
show the area covered by the FIR Bright Source Catalogue galax- 
ies. The redshifts of the sample are small and do not affect the 
colours. We observe that the FIR colour-colour relation in M 81 is 
within the consistent regime covered by the galaxy sample. Note 
that the FIR colours of the FIS Bright Source Catalog sample 
present the global colours, not those in individual regions within 
a galaxy. Thus, we confirm that the FIR colours of individual re- 
gions within a galaxy is fundamental in determining the global 
galaxy colours. The larger scatter of the FIS Bright Source Cata- 
logue sample may be due to a larger extent of the radiation field or 
a peculiarity of dust emission properties in some galaxies. 

In Fig. [6] we also show the colours of some representative re- 
gions within the circles of 1-arcmin radius as shown in Fig. [I] the 
central region ('C'), the spiral arms ('SI', 'S2', and 'S3'), and the 
interarm regions ('II' and '12'). We calculate the flux integrated for 
the circular regions, and take the flux ratios to show the colours. 
The fluxes are listed in TableQ] The central region have the bluest 
colours, while the interarm regions tend to have redder colours than 
the spiral arms: these trends in colours are consistent with the tem- 
perature map shown in Fig. [5] Moreover, the central and interarm 
regions are located relatively near to the theoretical predictions with 
varying ISRFs (squares in Fig. [6j on the colour-colour diagram, 
while the spiral arms are shifted toward the theoretical predictions 
with a mixture of general ISRF and a high radiation field (asterisks 
in Fig. [6j. This indicates that the ISRF in the central and interarm 
regions are rather uniform and occupied with general ISRFs with 
different intensities, while there is a large variation in the radiation 
field in the spiral arms with a spatial scale well below 1 arcmin. 
It is probable that small star-forming regions with strong radiation 
fields reside in spiral arms as seen in the Ha image (Fig.[T); this 
may be the reason for the mixed colour features of the spiral arms. 
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Figure 6. The relation between 65 ixm— 90 p.m colour and 140 ixm— 90 p.m 
colour. Small filled diamonds show the colours in the individual grids in 
M 8 1 . Triangles represents the colours for the flux integrated in the regions 
depicted in Fig. [T] Large open diamonds present the theoretical colour- 
colour relation with % = 1,3,10, 30, and 100, and asterisks show the results 
with / h = 0, 10~ 2 ' 5 , 10" 2 , 10" 1 ' 5 , 10 _1 , lO -0 ' 5 , and 1 from large to 
small 7"„(140 p.m)//, y (90 p.m). The area shown by the large square with 
dashed lines show the representative colour-colour relations of the galaxy 
sample in the AKARI FIS Bright Source Catalogue JPollo et alj201ol) . 

5 DISCUSSION 

5.1 FIR as a tracer of general ISRF 

The dust temperature derived from the global SED (Section l3.2t 
is consistent with those in the interarm regions and in the smooth 
regions (excluding the warm spots) of the spiral arms. This suggests 
that the global FIR SED is more representative of the dust heated 
by the general ISRF rather than by the intense radiation in the star- 
forming regions. 

That the FIR emission traces the general ISRF does not nec- 
essarily mean that the FIR traces the smoothly distributed old stel- 
lar population rather than the newly formed stellar population. It 
is gene rally hard to distinguish t hese t wo components as pointed 
out by IWalterbos & Greenawala n 199a) : An elaborated treatment 
of multi-wavelength sta r -form ation indicators can be found in 
Hiras hita. Buat. & Inouel J2003l) . who find that ~ 40 per cent of 
the FIR emission on average is associated with the dust heated by 
stellar populations older than 10 8 yr in their nearby star-forming 
galaxy sample. 

5.2 Colour-colour diagram as a diagnostics for the local 
heating source 

The dust optical depth derived observationally depends on the 
structures smaller than the spatial resolution. If the ISRF is uniform 
in a grid used to derive the optical depth, the dust optical depth is 
reliable because of the dust temperature estimated represents the 
real dust temperature. However, if a grid is contaminated by small- 
scale warm spots such as a star-forming clouds, the temperature 
estimate is biased to such warm spots. Because of this overestimate 
of dust temperature, the dust optical depth can be underestimated as 
pointed out in Section |4~T| In other words, the FIR emission is not 
a good tracer of dust mass if the region focused on contains small- 
scale warm regions contributing significantly to the FIR emission. 



We have shown in Section l4~2l that the contamination of warm 
regions can be identified by the FIR colour-colour diagram. As 
we can clearly see in Fig. [6] the spiral arms (SI, S2, and S3), 
which contain small-scale star-forming regions (see the Ha im- 
age in Fig. [TJ, deviate from the the varying-x sequence. The FIR 
colours in the spiral arms are more consistent with a mixture of 
a general ISRF with warm regions. Thus, we suggest that the rel- 
atively small dust optical depth in high-temperature regions of the 
spiral arms is due to the overestimate of dust temperature due to the 
contamination of warm spots possibly associated with star-forming 
regions (Section r4.lt . 

In the central region marked by 'C, the FIR colours are close 
to the point of uniform radiation field with \ ~ 10- Thus, the region 
may be roughly approximated with a single dust temperature. If this 
is true, the optical depth in the central region may reflect the real 
optical depth; that is, the deficiency of dust optical depth in the 
central region with a radius of 1 arcmin (1.1 kpc) may be real. The 
uniformity of the radiation field also implies that the heating by 
AGN is not a dominant heating source in the central 1-kpc region 
at least for the dust which can be traced i n the FIR bands. The H I 
gas is also deficient in the central region jRot j| 19751 : [Braunll 1995T) . 
Since the star formation activity in the central region is strong as 
seen in the Ha image, the radiation pressure or the thermal pressure 
may have pushed the gas and dust outwards. If we adopt 10 km s~ 
(~ the sound velocity in the ionized regions) for the typical velocity 
for the gas, a structure with a size of 1.1 kpc can be formed on a 
time-scale of 10 8 yr. Or some wave motion m ay be responsible for 
such a structure jRotsll 19751 : lLowe et al.lll994» . 



5.3 Importance of AKARI FIS bands toward Herschel 

The FIR wavelengths observed by the AKARI FIS bands (50- 
180 \im) cover the intensity peak of the FIR SED. This means that 
the intensities in the AKARI FIS bands depends both on the wave- 
length and the dust temperature in a strongly nonlinear way. This 
strong nonlinearity is important in the behaviour in the colour- 
colour diagram: not only the strong dependence on \ but also the 
large separation of the two sequences (i.e. varying \ an d varying 
/h) is due to such a nonlinearity. Moreover, the three (in fact four) 
bands available for AKARI FIS have advantage in the colour-colour 
analysis. 

Herschel is also suitable for the aim of mapping dust emis- 
sion in galaxies. In spite of its higher spatial resolution, the warm 
spots associated with star-forming regions cannot be resolved at 
the distance of M 8 1 . Thus, the only way to show the existence of 
such war m spots in the FIR i s to plot the FIR colours (Section l572l . 
Recently, iBendo et alj feOld) have observed M 8 1 with Herschel 
over a wide wavelength range of 70-500 \im. The peak of the SED 
is at A < 160 \im for T LG > 20 K. Thus, the Herschel PACS 
100 \im and 160 \xm bands (Poglitsch et alj|2010l) are suitable for 
covering the peak. The fact that similar temperatures are derived 
by both AKARI and Herschel (2 9 K, 20 K, 17 K f or the nucleus, 
arm, and interarm, respectively; IBendo et alj|2010h indicates that 
the FIR SED can be interpreted consistently from 90 [im to 500 
[im with a single temperature. Note that even with the Herschel 
higher resolution data, AKARI data still have an advantage that an 
enormous number of the AKARI All-Sky Survey data are available 
for a comparison at the same wavelengths (Section ^, 2l l. 
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6 CONCLUSION 

We have investigated the properties of FIR emission in a spiral 
galaxy, M81, by utilizing the AKARI imaging data at 65, 90, and 
140 [im. Combining the images in the two long-wavelength bands 
(90 and 140 \an), we have derived the dust temperature map. The 
dust temperature is ~ 25 K in the centre and becomes lower in 
the outer part. The dust temperature derived from the global 90 
and 140 vim intensities is Tlg = 18.6l}'g K, which reflects the 
dust temperatures in the interarm regions or the spiral arms exclud- 
ing the bright knots, rather than those in the centre or in the bright 
knots. Thus, the global dust temperature is more representative of 
the dust heated by the general ISRF. We have also shown the 140 
\im emission is more radially extended than the 90 \im emission, 
which is consistent with the radial dust temperature gradient. The 
dust optical depth traced by the FIR emission is Av = 0.5-2. If 
these values are divided by 2 with an assumption that the stars are 
in the mid-plane on average, the expected extinction for the stellar 
light is Av = 0.25-1, consistent with the the extinction derived 
from the FIR-to-UV ratio (Ay = 0.57). 

We have also demonstrated that the FIR colour-colour dia- 
gram is a useful tool to distinguish whether or not individual re- 
gions within a galaxy is dominated by a smooth ISRF or contam- 
inated by warm spots associated with star-forming regions. Based 
on this 'tool', we conclude that the bright regions in the spiral arms 
contain small-scale warm regions possibly hosting the H II regions 
seen in the Ha map. This contamination of warm regions causes an 
underestimate of dust optical depth (or dust column density). 

Since even Herschel cannot resolve individual star-forming re- 
gions, the FIR colour-colour diagram continues to be a useful tool 
to see if the dust is predominantly heated by a general ISRF or 
the dust heating by individual star-forming regions affects the FIR 
emission. 
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